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ABSTRACT 

We present an observational analysis of the numerical simulations of galaxy cluster 
mergers presented in Roettiger, Loken & Burns (1996). We identify several observational 
signatures of recent merger activity, and quantitatively assess the uncertainty introduced 
into cluster mass estimates when invoking the commonly held assumptions of hydrostatic 
equilibrium, virial equilibrium, spherical symmetry and isothermality. We find that 
mergers result in multiple X-ray peaks, long-lived elongation of the X-ray emission as well 
as isophotal twisting and centroid shifting to a degree consistent with recent observations. 
We also find an enlargement of the X-ray core relative to the dark matter core. Mergers 
result in non-isothermal clusters exhibiting observable inhomogeneities in the emission- 
weighted X-ray temperature of several keV on linear scales of less than 0.5 Mpc. The 
resulting gas dynamics are extremely complex, and we present an example of what might 
be observed by a high resolution X-ray spectrograph. We further speculate that the gas 
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dynamics, via shocks, bulk flows and turbulence, play an important role in the evolution 
of cluster galaxies and associated radio sources, particularly wide-angle tailed (WAT) 
sources and radio halos. We find that X-ray based cluster mass estimates made under 
equilibrium assumptions can be uncertain by 50% or more in the first 2 Gyrs after a 
merger and by up to 25% after 2 Gyrs depending on the details of the analysis and 
projection effects. Uncertainties can be considerably larger if the temperature is not 
well constrained. Similar uncertainties are observed in the X-ray derived baryon mass 
fractions. Virial mass estimates are typically overestimated because the observed 1- 
dimensional velocity dispersion can be severely contaminated by the infall velocity of 
the subcluster. 

Subject headings: Galaxies: clusters: general - galaxies: intergalactic medium - hydro- 
dynamics - methods: numerical - X-rays: galaxies 
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1. INTRODUCTION 

Clusters of galaxies are the largest gravitationally 
bound objects in the Universe. To date, thousands 
have been identified through both optical and X-ray 
surveys. Since they are abundant, luminous, and ex- 
tended objects capable of being observed in some de- 
tail at cosmologically significant distances, their value 
to cosmology is unquestioned. Clusters of galaxies 
have been used to constrain the matter density in the 
Universe, f2 (e.g., Nakamura, Hattori, & Mineshigc 
1995), the ratio of baryonic to dark matter density 
in the Universe, fif, (e.g., White & Fabian 1995), and 
even the power spectrum of primordial density fluc- 
tuations (e.g., Henry & Arnaud 1991). In each case, 
accurate knowledge of the gravitational mass in clus- 
ters is required. Traditionally, cluster masses have 
been estimated from the velocity dispersion and spa- 
tial distribution of member galaxies assuming them 
to be in virial equilibrium (e.g., Merrit 1987), or from 
the distribution of X-ray emitting gas assuming it to 
be in hydrostatic equilibrium (e.g., Fabricant, Lecar, 
& Gorentcin 1980). More recently, it has become pos- 
sible to determine cluster masses from the effects of 
gravitational lensing (e.g., Tyson, Valdez, & Wenk 
1990). Since this is the only direct method of mass 
determination and is therefore not subject to assump- 
tions regarding cluster dynamics, it is potentially the 
most powerful tool. Unfortunately, in the instances 
where all three methods have been applied, the results 
can often disagree by factors of several thus raising 
questions as to the validity of the underlying assump- 
tions in each method (e.g., Loeb & Mao 1994; Kncib 
ct al. 1995). 

In a hierarchical universe in which large-scale struc- 
tures form via the accretion of smaller structures, 
the greatest threat to equilibrium conditions within 
a given cluster resides in the potential for a merger 
with another cluster. In order to determine the degree 
and duration of non-equilibrium conditions following 
a major merger event, we have conducted a series 
of numerical hydrodynamical/N-body simulations of 
isolated merging clusters of galaxies. The results of 
these simulations have appeared in a series of papers. 
Paper I (Roettiger, Burns & Loken 1993) presented 
a brief description of the numerical simulations along 
with several general results while Roettiger, Loken & 
Burns 1996 (Paper II) provided a detailed description 
of the numerical simulations and the physical evo- 
lution of clusters during the merger. In this paper, 



we address the questions: 1) What are the observ- 
able signatures of merging clusters? 2) How does the 
violation of the conditions of hydrostatic and virial 
equilibrium affect the determination of such cosmo- 
logically important parameters as cluster masses and 
the baryon to dark matter mass ratio? Understand- 
ing the observational signatures of cluster-subcluster 
mergers will help us to identify clusters which are not 
relaxed and may therefore be subject to considerable 
uncertainty in their mass estimates. 

The numerical simulations and the parameter sur- 
vey upon which this work is based are discussed in 
detail in Paper II. Here, we present only a brief 
overview. The numerical simulations employ a hy- 
brid Hydro/N-body code in which the intracluster 
medium is represented by the equations of hydrody- 
namics, and the dark matter is represented by the N- 
body particle distribution. The hydrodynamical com- 
ponent of the code is ZEUS-3D (Clarke 1990; Stone 
& Norman 1992), an Eulcrian, finite-difference code 
developed and supported by the Laboratory for Com- 
putational Astrophysics (LCA) at the National Cen- 
ter for Supercomputing Applications (NCSA). The N- 
body component of the code is TREECODE (Hcrn- 
quist 1987), a hierarchical tree algorithm in which the 
force between nearby particles is calculated directly 
while more distant particles are grouped so that the 
force term is calculated via a multipole expansion. 
The gravitational potential is defined entirely by the 
dark matter component. The gas component is cho- 
sen to be <5% of the total cluster mass so as to be 
dynamically insignificant. Radiative cooling has not 
been included in these simulations. The central clus- 
ter densities, nominally 10 -3 cm~ 3 , were chosen so as 
to make the effects of cooling negligible. 

We begin with two isolated, isothermal, hydro- 
static King spheres placed on the computational grid 
such that their cores are separated by ~5 Mpc. They 
are then allowed to merge under the influence of their 
mutual gravity. The idealized initial conditions have 
several benefits. First, they allow a more efficient 
use of the computational volume, and consequently, 
higher resolution than is typically found in most cos- 
mological simulations (e.g., Bryan et al. 1994). The 
resolution of these simulations is 50 kpc or 5 zones 
across the primary cluster core radius. Second, they 
provide a well-defined baseline with which to com- 
pare the subsequent evolution of the clusters after 
the mergers. And third, they allow the study of the 
merger event in the absence of the more general clus- 
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ter formation process. 

The parameter survey discussed in Paper II was 
designed to study the effects of varying the relative 
cluster-subcluster masses and their relative gas con- 
tent. The mass ratios in the survey were 1:1, 2:1, 
4:1 and 8:1. The variation of the cluster-subcluster 
mass ratio implied different values for both the im- 
pact velocity and the relative linear dimensions of the 
clusters. We also conducted a simulation in which 
the smaller of the two clusters (hereafter, the sub- 
cluster) was essentially gas-free. This showed the 
relative importance of the hydrodynamical processes 
(i.e., shocks, bulk flows) in the merger evolution. 

In this paper, we address some of the observational 
consequences of cluster-subcluster mergers. In {S, 
we describe the merger samples presented in this pa- 
per. In §^], we discuss the evolution of the X-ray 
morphology. Section |I] contains a discussion of the 
observable temperature substructure. Section [[] ad- 
dresses the observable consequences of the gas dynam- 
ics. Section [3] examines the errors associated with 
cluster mass determinations based on the commonly 
employed assumptions of spherical symmetry, hydro- 
static and virial equilibrium. Also in ^6], we discuss 
the implications of uncertain mass estimates and non- 
equilibrium conditions on the determination of the 
baryon fraction in clusters of galaxies. Section |7] is a 
summary of our results. 

2. THE SAMPLES 

For the sake of brevity and clarity, we cannot 
present every observational aspect of every merger in 
our parameter survey. We have therefore chosen two 
sub-samples having characteristics representative of 
the whole. The first is the evolutionary sample which 
shows the evolution of the 4:1 mass ratio merger at 
six epochs. The second, or deprojection sample, con- 
tains several mergers representing various mass ratios, 
projection angles and epochs. The chosen epochs are 
-0.25, 0, 0.5, 1.0, 2.0 and 5.0 Gyrs relative to the time 
of core passage. These epochs are the same as those 
presented for all mergers in Paper II thus allowing a 
direct examination of the underlying physics at each 
epoch. 

2.1. The Evolutionary Sample 

The evolutionary sample was chosen to depict the 
evolution of a single mass ratio merger as viewed in 
the plane of the sky. We have chosen the 4:1 mass ra- 



tio merger as a representative example of a moderate 
mass ratio merger. At each epoch, we view it perpen- 
dicular to the merger axis because it is from this angle 
that the X-ray morphology is most severely disturbed. 
Observational projection effects are addressed in the 
deprojection sample (§2.2.). To summarize, the pri- 
mary cluster is scaled to 4 x 10 14 M Q with a velocity 
dispersion of 660 km s _1 . The sub-cluster is scaled to 
1 x 10 14 M© with a velocity dispersion of 415 km 
s . The primary cluster gas temperature is 4.6 keV 
while the subcluster temperature is 1.8 keV. The dark 
matter component of the clusters is represented by 
22,500 equal mass particles. The synthetic X-ray im- 
ages for this sample can be seen in Figure (See §||]). 

2.2. The Deprojection Sample 

The deprojection sample was chosen to be repre- 
sentative of various mass ratios, epochs, and merger 
axis projection angles such that a variety of realis- 
tic and representative observational situations are de- 
picted. To this end, we have selected six examples for 
the X-ray deprojection analysis (§^|). The synthetic 
X-ray images can be seen in Figure 0. The merger 
parameters as well as several observable quantities 
are listed in Table [|. Further details are in Table 
1 of Paper II. The sample includes two examples with 
obvious X-ray substructure (Mergers 1 and 5, as in- 
dicated in Table 1), three examples in which the X- 
ray morphology appears somewhat relaxed although 
it is known that shocks are present (Mergers 2,3, and 
4), several mergers in which the observed velocity 
dispersion is contaminated by the infall velocity to 
varying degrees (Mergers 1-4), and two examples in 
which shocks and bulk gas flows have diminished and 
contamination of the velocity dispersion is minimal 
(Mergers 5 and 6). 

3. THE X-RAY MORPHOLOGY 

In this section, we examine the variety and evo- 
lution of the X-ray surface brightness morphologies 
exhibited during the merger. Synthetic X-ray im- 
ages of the evolutionary and deprojection samples can 
be seen in Figures [j] and |[ respectively. The X-ray 
emission is based on the Raymond-Smith X-ray spec- 
trum which includes both thermal bremsstrahlung 
and line emission within a ROSAT-like bandpass of 
0.5-2.0 keV (Raymond & Smith 1977). Using the 3- 
dimensional density and temperature data from the 
numerical simulations, we have calculated the X-ray 
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volume emissivity in each zone. The X-ray surface 
brightness images (Figures and 0) were generated 
via a line-of-sight (LOS) integration through the com- 
putational volume assuming the emission to be opti- 
cally thin. No attempt has been made to add statis- 
tical noise. The boxy edges in the outer contours are 
indicative of the reduced resolution near the edges of 
the computational volume (see Paper II, §3.4). 

Figure shows the evolution of the X-ray surface 
brightness for the evolutionary sample. Not pictured 
is the initial cluster which was spherically symmetric. 
The surface brightness has been projected such that 
the merger axis runs vertically and is in the plane of 
the sky. The subcluster entered from the top of the 
panel and proceeded downward. At the top of each 
panel, we indicate the merger epoch relative to the 
time of core passage. The contour levels are logarith- 
mic and constant from one epoch to the next. 

In Figure 0, we see a variety of X-ray morpholo- 
gies. Several epochs, both before and after core pas- 
sage, exhibit multiple concentrations of X-ray emis- 
sion (Figures 0a,c,d,e). In this merger, the subcluster 
penetrated the core of the primary and exits on the 
downstream side resulting in a secondary maximum 
below the X-ray core. All epochs after core passage 
exhibit a general elongation of the X-ray morphology 
parallel to the merger axis. As was shown in §5.5 of 
Paper II, the elongation of the X-ray (i.e., gas density) 
distribution is often more extreme than that of the un- 
derlying dark matter distribution. We also note that 
the elongation of the X-ray morphology is not limited 
to being parallel to the merger axis. Figures 0b and d 
and to a lesser extent, e, show distinct elongation of 
the central region perpendicular to the merger axis. 
Although not apparent in Figure 0c because of the 
choice of contours, both epochs corresponding to Fig- 
ures 0c and d show "bars" of X-ray emission near the 
cluster core. An equivalent feature is not apparent in 
the dark matter distribution indicating a severe vio- 
lation of hydrostatic equilibrium. The elongation of 
the gas distribution perpendicular to the merger axis 
is due to a shock and adiabatic compression gener- 
ated by the infalling subcluster. This morphology is 
quite similar to that of A754 (Henry & Briel 1995), 
and Henriksen (1993) has suggested that A754 is the 
result of a recent merger between two unequal mass 
clusters. 

Finally, even after 5 Gyrs (Figure 0f), we still see 
the effects of the merger in the X-ray morphology. 
On the downstream side of the core there is a flatten- 



ing of the X-ray contours as a result of residual infall 
from the subcluster 's return trip through the primary 
core. Also notable in Figure 0f is the much expanded 
X-ray distribution (compare with Figure 0a). One 
important consequence of the merger is an increase 
in the cluster X-ray core radius, often by a factor of 
2, relative to the dark matter core. The expansion of 
the gas core is largely thermally supported as a re- 
sult of shock heating. However, in the early stages 
of the merger dynamical support (i.e., bulk flows and 
turbulence) is also important. At 0.5 and 1.0 Gyr, 
the mean dynamical pressure (pv 2 ) within the core 
is nearly 30% of the mean thermal pressure (nkT) 
within the same region. At 2.0 and 5.0 Gyrs the mean 
dynamical pressure has dropped to less than ~15% of 
the mean thermal pressure. Figure demonstrates 
both the absolute expansion of the X-ray core as well 
as the expansion of the X-ray core relative to the dark 
matter core. 

The significance of the relative core sizes is two- 
fold. First, recent gravitational lensing experiments 
find that dark matter cores are significantly smaller 
than the X-ray cores (e.g., Miralda-Escude & Babul 
1995) leading to lensing-based mass estimates 2-3 
times larger than those derived from the X-ray emis- 
sion. Second, the standard X-ray deprojection tech- 
nique (Arnaud 1988) requires an estimate of the 
dark matter or gravitational core radius. Since the 
dark matter is unobservable, the core radius estimate 
comes from the X-ray surface brightness image. The 
consequences of this discrepancy are discussed in § |fi] . 

Figure contains the synthetic X-ray surface bright- 
ness images for the deprojection sample. Again, the 
epochs depicted are the same as in Figure [l] only now 
the mass ratios and merger axis projection angles dif- 
fer. The X-ray morphologies shown in Figure are 
quite similar to those in Figure 0. However, the pro- 
jection can often serve to soften features and hide 
evidence of a recent merger. This is important since 
the X-ray emission is often used as an indicator of 
the evolutionary state of the cluster (e.g., Mohr et 
al. 1995). 

One of the best examples of how projection effects 
can hide substructure is Figure 0d in which a 2:1 mass 
ratio merger is projected 45° from the plane of the sky 
1 Gyr after core passage. If viewed in the plane of the 
sky, one would witness a considerably more prominent 
extension below the cluster core as well as extensions 
perpendicular to the merger axis on either side of the 
core (e.g., Figure |l]i). These features are almost com- 
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pletely hidden in this projection. In order to obtain 
perspective on the projection effects, one can compare 
the projected X-ray morphologies in Figure^] with the 
gas density contours for the corresponding mass ratio 
and epoch in Figures 6-9 of Paper II. 

In the absence of obvious substructure, we can look 
for more subtle evidence of merger activity such as 
isophotal twisting or centroid shifting. The results of 
this test are shown in Figure ^ where we have used the 
IRAF routine ELLIPSE to fit elliptical isophotes to 
the synthetic X-ray distribution in Figure ^d. In this 
case, the ellipticity varies dramatically with radius 
from near at large radii to 0.42 near the core. Al- 
though there does appear to be some minor isophotal 
twisting, all position angles of all ellipses are consis- 
tent at the 2-cr level with alignment along the merger 
axis. We also see statistically significant centroid 
shifts of ~150 kpc. Nearly the entire magnitude of 
the centroid shift occurs along the merger axis. The 
magnitude of the maximum centroid shift in this ex- 
ample is only slightly larger than the largest reported 
by Mohr et al. (1995) for their sample of clusters 
observed by the Einstein satellite. 

4. THE GAS TEMPERATURE MORPHOL- 
OGY 

The gas temperatures within the ICM can be an 
important diagnostic of the dynamical state of the 
cluster. Substructure, shocks and adiabatic cool- 
ing can result in a very patchy and distinctly non- 
isothermal gas distribution. This is in stark contrast 
to the commonly held notion of an isothermal cluster. 
Because of limited spatial and spectral resolution in 
previous X-ray observations, it is only recently that 
the level of inhomogeneity in the temperature mor- 
phology has become fully appreciated (e.g., Briel & 
Henry 1994; Arnaud et al. 1994). The situation, how- 
ever, is rapidly changing since ASCA, and, eventually 
AXAF, promise to supply even more accurate and 
detailed temperature maps in the near future. For 
this reason, we focus here on the observable tempera- 
ture morphology that results from cluster-subcluster 
mergers. 

Figures || and || show the observable temperature 
maps (solid contours, labeled with the temperature in 
keV) overlaid onto the X-ray surface brightness (dot- 
ted cont ours ) for the evolutionary (§2.1.) and depro- 
jection §2.2.) samples, respectively. The temperature 
maps are generated by calculating a mean emissivity- 



weighted temperature along the LOS through the 
cluster. The temperature in each zone is weighted by 
the volume emissivity of that zone within a ROSAT- 
like bandpass of 0.5-2.0 keV. The true cluster gas tem- 
peratures can be seen in Figure 12-15 of Paper II. 
Comparing these figures shows that projection effects 
can severely alter the temperature morphology. In 
general, the projection and LOS integration tend to 
mitigate the extreme temperatures seen in the various 
mergers. 

As an example of projection effects, Figure |^a 
shows the 2:1 mass ratio merger shortly before core 
passage and projected at 45° to the plane of the 
sky. This is a similar situation to the A2256 model 
presented in Rocttiger, Burns, & Pinkney (1995) al- 
though the scaling is different. The true temperature 
in a central slice along the merger axis is shown in 
Figure 0. Close examination shows the relatively cool 
subcluster (top) impinging on the somewhat warmer 
primary cluster (bottom). A shock has formed at 
the cluster-subcluster interface and is visible as a thin 
layer of very hot gas perpendicular to the merger axis. 
This region is actually a broad conical structure whose 
axis of symmetry coincides with the merger axis. The 
hottest gas is located at the apex of the cone. In Fig- 
ure ^|a, the conical region of hot gas appears as two 
"hot lobes" located on either side of the merger axis 
(vertical line through the X-ray maxima) . Because of 
the 45° projection, the conical structure is seen more 
head-on and is therefore broader than the more linear 
feature seen in Figure [?]. It takes on the lobe structure 
because the temperature in the central region is con- 
taminated by cooler subcluster gas projected along 
the same LOS. Figure exhibits the expected mor- 
phology when the conical shock front is viewed in the 
plane of the sky. 

A variety of complex temperature morphologies 
are exhibited throughout the rest of the deprojection 
sample. In Figure Ifjjs, there are four distinct temper- 
ature regions located along the merger axis. Starting 
at the top and moving downward, there is a region 
contained within 7.5 keV. This region is actually a 
local minimum (<7.5 keV) and is identified as the 
relatively cool wake of the infalling subcluster. There 
is a hotspot (>16 keV) located at the X-ray core. Be- 
low the X-ray core, there is another cool region that 
we identify as the emerging subcluster. Below the 
subcluster there is a large arc of shock-heated gas, > 
11 keV. A similar morphology is seen in Figure ||d 
though in this case, the hotspot, once located near 
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the core, is seen to have propagated upstream. This 
is the result of a shock that forms when gas that is 
expelled from the core by the rapidly varying gravi- 
tational potential impinges on residual infall from the 
subcluster. To summarize, we find regions of gas vary- 
ing by several keV that are separated on linear scales 
of several hundred kiloparsecs. 

As one might expect, the most complex tempera- 
ture structures will be those viewed perpendicular to 
the merger axis and within 2 Gyrs after core passage. 
Figure ^e reveals a hot core which, to varying degrees, 
results from shock heating and a continuous stirring 
of the gas core by the varying gravitational potential. 
It also shows a significantly cooler region downstream 
from the merger. This appears to be analogous to a 
cool region noted by Schindler (1996) which was iden- 
tified as a rarefaction wave following the initial bow 
shock. Although this is a possible explanation, there 
is evidence in our simulations that this is actually gas 
from the subcluster that has passed through the core. 

As described in §5.2.1 in Paper II, we believe that 
the subcluster gas is protected from ram pressure 
stripping by the bow shock that forms on the subclus- 
ter's leading edge shortly before core passage. The 
effect of the bow shock is to reduce the ram pres- 
sure by reducing and redirecting the flow relative to 
the subcluster core. We tested this interpretation by 
running the 4:1 mass ratio simulation with a passive- 
scalar, (a dynamically inert quantity used to trace the 
various gas components in numerical simulations, see 
§5.2, Paper II), attached to the subcluster and found 
that the passive-scalar did penetrate the core. 

Since the gravitational potential of the subcluster 
is largely disrupted after core passage, it no longer 
confines the subcluster gas, and consequently, the 
subcluster gas expands and cools adiabatically. 

Finally, Figure |^f is viewed perpendicular to the 
merger axis but at this stage, 5 Gyrs after core pas- 
sage, we find a relatively relaxed configuration in 
which the shocks have dissipated and the subcluster 
is thoroughly mixed. The centrally peaked tempera- 
ture distribution is maintained by a continual mixing 
of the gas by the varying gravitational potential min- 
imum. 

The results presented here are quantitatively very 
similar to recent temperature maps generated from 
both ROSAT and ASCA data. We have made a 
detailed comparison between these models and the 
i?05'j4Ttemperature map of A2256 presented by Briel 



& Henry (1994) in Roettiger et al. (1995). Henry 
& Briel (1995) present a temperature map of A754 
which shows inhomogeneities of similar scale. Ini- 
tial ASCA results have been published for A2163 
(Markevitch et al. 1994) and Perseus (Arnaud et 
al. 1994). Both show widely varying temperatures 
(AT ~ 5 keV) on scales <0.5 Mpc. In addition to 
the temperature morphology, there is other evidence 
that A2163 is in the process of a merger. Miralda- 
Escude & Babul (1995) find evidence from gravita- 
tional lcnsing that the gas distribution is not in hy- 
drostatic equilibrium with the gravitational potential, 
and Herbig & Birkenshaw (1994) have identified a ra- 
dio halo associated with A2163. Cluster radio halos 
are a Mpc-scale phenomenon which has been linked to 
cluster-subcluster mer gers (J affe 1992; Tribble 1993; 
Burns et al. 1994, see |5.2.2.[). 



5. THE GAS DYNAMICS 

There are several important observational conse- 
quences of the gas dynamics induced by the cluster- 
subcluster merger. As we have seen, mergers result in 
shocks, bulk flows, and large eddies which appear to 
decay into a random velocity field on timescales of ~2 
Gyrs. The indirect effects of the gas dynamics on the 
accuracy of equilibrium-based mass estimates will be 
addressed in §p] . In this section, we discuss the more 
directly observable consequences. 

5.1. Spectroscopic Gas Velocities 

Currently, it is not possible to directly measure the 
dynamics of the X-ray emitting ICM. However, this 
will change with the planned launch of ASTRO-E in 
2000. So, with the launch date rapidly approaching, 
it is now time to make predictions as to the observable 
velocity field that may result from a cluster-subcluster 
merger. In this section, we use the results from the 
numerical simulation to determine quantitatively the 
level and structure of the observable gas dynamics 
as they might be viewed by a high spectral/spatial 
resolution X-ray telescope. 

As with other astronomical spectroscopy, X-ray 
spectroscopy will only allow the determination of LOS 
velocities and, assuming optically-thin emission, the 
X-ray spectrum at any point will be a composite of 
all emitting regions along the LOS. We have not, at 
this time, attempted a detailed model of the clus- 
ter X-ray spectrum nor have we taken into account 
the specific characteristics of ASTRO-E. We have, 
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however, attempted to model an "observable" veloc- 
ity map based on the projected LOS gas velocities 
weighted by the corresponding Raymond-Smith X- 
ray emissivity. Since the X-ray spectral lines are re- 
combination lines, the line emissivity, like that of the 
thermal bremsstrahlung, is expected to scale as n 2 
(e.g., Sarazin 1986). 

As an example of an observable velocity map, Fig- 
ure ||a is a subregion of Figure ||cl and contains the 
mean LOS emission-weighted gas velocity for Merger 
4 of the deprojection sample (see Table |l|). Similarly, 
Figure |^b contains the emission- weighted gas velocity 
dispersion along the same projected LOS. All veloci- 
ties are in the rest frame of the dark matter core of 
the primary cluster. As a spatial reference, we have 
included in Figure || contours of the X-ray surface 
brightness (dotted line). 

Merger 4 is a 2:1 mass ratio merger at 1 Gyr after 
core passage. The merger axis is projected 45° out of 
the plane of the sky such that the subcluster is falling 
toward the observer (negative velocity). The region of 
relatively large negative velocity, < —450 km s _1 , lo- 
cated below the X-ray peak is subcluster gas that has 
passed through the core of the primary. As with all 
other figures, the subcluster originated at the top of 
the panel and proceeds downward. The effect of the 
merger axis projection and the LOS integration is to 
soften the velocity extremes. The true velocity struc- 
ture, in a plane parallel to the merger axis, can be 
seen in Figure 25d of Paper II. The peak velocity in 
that figure is in excess of 1000 km s _1 . In Figure ||a, 
the absolute maximum is ^550 km s _1 . The effect 
of the emission-weighting is to bias the LOS integra- 
tion toward the denser and therefore more dynami- 
cally (and observationally) significant central regions 
of the cluster. 

The velocity dispersion map is also emissivity- 
weighted. We find the greatest velocity dispersions 
to be near the edges of the advancing subcluster (two 
peaks at 350 km s" 1 below the X-ray core) where 
gas is being sheared away and rapidly decelerated. A 
large dispersion is also apparent in the region above 
the X-ray core where there is still some residual infall 
that is being decelerated as it collides with gas ex- 
pelled from the core. It should be kept in mind that 
this is only one example and is subject to strong pro- 
jection effects. It is, therefore, meant only to serve as 
an illustrative example. 



5.2. Gas Dynamics and Cluster Radio Sources 

5.2.1. The WAT Problem 

Generally, radio tails are found to be underpressur- 
ized with respect to the ICM and are therefore sub- 
ject to the dynamical influences of the thermal X-ray 
gas (Feretti, Perola, & Fanti 1992). It has long been 
suggested that the narrow-angle tailed (NAT) radio 
source morphology results when the radio tails are 
swept back by ram pressure as the host galaxy moves 
at high speed through the ICM (Miley et al. 1972). 
Although this explanation works well for NATs, it 
does not work well for wide-angle tailed radio sources 
or WATs (e.g., Eilek et al. 1984). WATs are C-shaped 
radio sources consisting of two, straight and narrow 
jets that flare into broad tails which, like NATs, ex- 
tend in a common direction. The jets typically ex- 
tend 10-100 kpc while the total extent of the source, 
including tails, can reach 1 Mpc. Unlike NATs, how- 
ever, WATs are associated with the central dominant 
galaxies in clusters and are therefore presumed to be 
at rest with respect to the cluster gravitational poten- 
tial and, assuming hydrostatic equilibrium, the ICM. 
What then is bending the WAT? (see also Eilek et 
al. 1984; O'Donoghue, Eilek & Owen 1993) 

We have suggested that bulk flows generated dur- 
ing mergers may be responsible for the bending of the 
WAT tails (Roettiger et al. 1993; Pinkney, Burns & 
Hill 1994; Loken et al. 1995). The numerical sim- 
ulations presented here allow us to add quantitative 
support to this idea. Previous analysis of the radio 
data indicates that ICM flow velocities on the order 
of 1000 km s" 1 are required, assuming typical central 
gas densities (~10~ 3 cm~ 3 ), to provide the necessary 
ram pressure to bend the jets (e.g., Burns et al. 1986). 
The linear scale of WATs, ranging from 100 kpc to 1 
Mpc, gives some indication as to the required linear 
dimensions of the bulk flows. These conditions are 
met for extended periods during the evolution of a 
merger. 

We find peak gas velocities well in excess of 1000 
km s _1 at various stages of the merger evolution. Fig- 
ure H shows that the peak gas velocities generally do 
not decay below 1000 km s _1 for nearly 2 Gyrs after 
core passage. This is much longer than the estimated 
ages of the radio sources, 10 7 -10 8 years. Figure 10 
shows the level of coherence to the flow as a function 
of time. Typically the bulk flow through the clus- 
ter core is >200 kpc wide and can be ~1 Mpc wide. 
The coherence of the flow diminishes on time scales of 
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several Gyrs. Finally, we find bulk flows in both di- 
rections along the merger axis at various times during 
the evolution of the merger. 

If the merger hypothesis is true, one would ex- 
pect WAT clusters to exhibit other merger signatures 
such as an elongation of the X-ray morphology in the 
same direction as the radio tails, an absence of cooling 
flows, velocity substructure, etc. All of these are in- 
deed observed. Velocity substructure in WAT clusters 
is discussed in Pinkney (1995). ROSAT observations 
showing elongations of X-ray emission along the di- 
rection of the WAT tails are presented in Gomez et 
al. (1996) and in Pinkney et al. (1993, 1994). Nu- 
merically, this hypothesis has been addressed more 
directly by Loken et al. (1995) with numerical sim- 
ulations of radio jets/tails within the velocity field of 
a merging cluster like that in Figure 24c of Paper II. 

5.2.2. Radio Halos 

Radio halos are large (~ 1 Mpc), amorphous clus- 
ter sources of uncertain origin and generally steep ra- 
dio spectra (see Hanisch 1982). They are extremely 
rare sources that appear to be associated only with 
very rich clusters that have undergone recent merg- 
ers. Some of the best examples include Coma, A2255 
and A2256. Each of these clusters has been reported 
to have undergone a recent merger based on a variety 
of properties including X-ray morphology, galaxy ve- 
locity distributions, temperature morphology and the 
absence of cooling cores (Burns et al. 1994, Burns 
et al. 1995, Roettiger et al. 1995, Briel et al. 1991). 
Therefore, as has been suggested by various authors 
(e.g., Tribble 1993; Burns et al. 1994), it may be rea- 
sonable to expect that cluster halos are also related 
to recent merger activity. 

The difficulty in explaining radio halos arises from 
the combination of their large size, up to 1 Mpc, and 
the short synchrotron lifetime of the relativistic elec- 
trons. The expected diffusion velocity of the electron 
population is on the order of the Alfven speed (~100 
km s _1 ) making it difficult for the electrons to dif- 
fuse over a Mpc-scale region within their radiative 
lifetime. Therefore, one needs a mechanism by which 
the relativistic electron population can be transported 
over large distances in a short time, or a mechanism 
by which the local electron population is reacceler- 
ated and local magnetic fields are amplified over an 
extended region. The cluster-subcluster merger can 
potentially supply both mechanisms. 



The first mechanism has been used to explain the 
halos in A2256 (Rottgering et al. 1994) and Coma 
(Giovannini et al. 1993). Rottgering et al. (1994) 
contend that the halo source in A2256 results from 
the distortion of several pre-existing head-tail sources 
by the turbulence and shear generated by the infalling 
subcluster. They calculate that the relativistic elec- 
trons would need to be transported at 1500-3000 km 
s _1 in order to account for the extent of the source. 
In the A2256 merger model presented in Roettiger et 
al. (1995), the infalling subcluster has a velocity of 
nearly 3000 km s _1 at the current epoch. 

The second mechanism has been applied to the 
Coma cluster radio halo (Burns et al. 1994). Here, 
the radio halo results from the global disruption 
caused by the merger. During the merger, shocks and 
turbulence create sites for in-situ particle reaccelera- 
tion and amplification of the magnetic field over large 
scales (Eilek & Hughes, 1990; Ruzmaikin et al. 1989). 
A simple calculation of the total kinetic power im- 
parted by the merger (e.g., Burns et al. 1995) shows 
that the merger can potentially supply many orders 
of magnitude more energy than is necessary to power 
the radio source, 10 46 vs. 10 41 ergs s _1 for Coma 
and A2255 (Burns et al. 1994; Burns et al. 1995). 
Although the details of this mechanism are not fully 
understood, if we accept this premise, the lifetime of 
the halos will be governed not by the radiative lifetime 
of the electron population but rather by the time scale 
for the dissipation of the shocks and decay of the tur- 
bulence. Looking at the bulk velocity properties pre- 
sented in Figures ^ and [LlJ we find that the e-folding 
time for both the peak and mean gas velocities asso- 
ciated with the central 1 Mpc of the primary cluster 
is greater than 1 Gyr. The shocks dissipate on sim- 
ilar timescales. This opens the possibility that radio 
halos may live considerably longer than the 10 8 yrs 
suggested by Tribble (1993). As Tribble noted, one 
benefit of the short lifetime is that it serves to explain 
the rarity of radio halos. However, their rarity may 
be attributable to the relative rarity of suitable hosts. 
To date radio halos have only been identified in very 
rich clusters (richness class >2). Finally, if both pro- 
posed mechanisms are responsible for the generation 
of halos, it may be that each is important during dif- 
ferent stages of the merger evolution. A2256 appears 
to be in the very early stages of its merger (Roettiger 
et al. 1995) while Coma is in the relatively late stages 
of its merger (Burns et al. 1994). 
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5.3. Gas Dynamics and Galaxy Evolution 

The galaxies found in clusters of galaxies are sys- 
tematically different from those in the field (see Oem- 
ler 1992 for a review). Cluster galaxies are typically 
redder and contain less gas than isolated field sys- 
tems. They are also more likely to be ellipticals than 
spirals, particularly in the cores of rich clusters. This 
has led many to argue that the cluster environment 
plays an important role in the evolution of individual 
galaxies (e.g., Spitzer & Baade 1951; Gunn & Gott 
1972). 

One mechanism by which the environment may in- 
fluence galaxy evolution is via the galaxy-ICM in- 
teraction. Dressier & Gunn (1983) suggested that 
field galaxies fall into the cluster environment where 
they encounter shocks in the ICM. As a result of 
the shock encounter, the galaxies experience an in- 
crease in ram pressure which may induce a burst of 
star formation and/or strip the galaxy of its inter- 
stellar medium (ISM) . Numerical simulations of clus- 
ter formation show that during the formation pro- 
cess a spherical accretion shock forms near the core of 
the cluster and propagates outward effectively form- 
ing an outer boundary to the cluster (Evrard 1990). 
Evrard (1991) showed that field galaxies undergoing 
a burst of star formation upon interacting with the 
accretion shock may explain many of the properties 
of Butcher-Oemler galaxies observed at moderate red- 
shifts (Butcher & Oemler 1978). Edge et al. (1990) 
suggested that cluster mergers may also contribute to 
the Butcher-Oemler effect since the apparent galaxy 
evolution occurs during a period of significant evo- 
lution in the cluster X-ray luminosity function (z < 
0.5). 

A bow shock forms on the leading edge of the 
subcluster as it approaches the core of the primary. 
This shock serves to protect the subcluster from be- 
ing stripped of its own ICM by creating a boundary 
layer that effectively reduces the ram pressure by de- 
creasing the flow velocity and redirecting it around 
the core of the subcluster (see Figure 10 of Paper II). 
In doing so, the shock will also serve to protect the 
relatively gas-rich subcluster galaxies. The protection 
fails, however, as the cores become coincident. As the 
shock approaches the core of the primary, it encoun- 
ters a rapidly increasing density profile and deceler- 
ates with respect to the subcluster dark matter and 
presumably the subcluster galaxy population. This 
allows the galaxies to pass through the shock initiat- 



ing a burst of star formation followed by a rapid strip- 
ping of their ISM which effectively truncates the star- 
burst. It is also possible that the subcluster galaxies 
are already actively forming stars and that interaction 
with the bow shock simply truncates the existing star 
formation. In either case, it is this sort of truncated 
star formation that is believed to be the mechanism 
by which post-starburst galaxies are formed (Dressier 
& Gunn 1983; Newberry, Boroson & Kirshner, 1990). 

Using our numerical simulations and assuming that 
the dark matter particle dynamics are representative 
of galaxy dynamics, we can track the ram pressure 
history of individual particles as a function of their 
initial cluster membership. We find that the subclus- 
ter particles, on average, do experience an impulsive 
burst of ram pressure as the cluster cores merge (see 
Figure [l2|) . In fact, the median ram pressure experi- 
enced by subcluster particles increases briefly by fully 
two orders of magnitude over the initial value. We 
have already used this mechanism to explain the post- 
starburst galaxies observed in the bridge between 
Coma and the NGC 4839 galaxy group (Burns et 
al. 1994; Caldwell et al. 1993). These galaxies are 
dispersed in a manner consistent with the subcluster 
N-body particles in the merger simulation. That is, 
they are distributed along an extension in the X-ray 
emission connecting Coma and the NGC 4839 group. 
For the higher mass ratio mergers, as we believe Coma 
has recently undergone, the induced ram pressure is 
very impulsive. We were therefore able to use the 
short duration of the burst in conjunction with star- 
burst evolution models to constrain the age of the 
merger. 

In more general terms, our simulations show that 
cluster-subcluster mergers can potentially introduce 
expansive shocks into the cluster environment. This 
is quite apparent from Figure O. In this Mach num- 
ber contour plot, the shocks are indicated by strong 
gradients or closely spaced contours. Generally, we 
see two shocks form as a result of the merger. They 
often appear as two opposing C-shaped structures in 
the Mach number figures. The first shock to form 
is the initial bow shock described above. This shock 
propagates downstream while expanding and eventu- 
ally dissipating. A second shock, as described in §[h], 
forms after core passage and propagates upstream. 
As is apparent in Figure these structures can be 
of Mpc-scale and can sweep through Mpc-scales be- 
fore completely dissipating and thus may influence 
the evolution of many subcluster and primary cluster 
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galaxies. 

6. EQUILIBRIUM MASS ESTIMATES 

The various methods used by observers to estimate 
the masses of galaxy clusters often require the as- 
sumptions of hydrostatic equilibrium, virial equilib- 
rium, and spherical symmetry. In this section, we at- 
tempt to reproduce observational conditions and test 
these methods on the simulated clusters during vari- 
ous stages of the merger evolution. Under equilibrium 
conditions, the total mass enclosed within a radius, r, 
is given by Fabricant et al. (1980) as, 



M(r) 



kT{r)r 
G\im v 



dlnp(r) dlnTir) 



dlnr 



dlnr 



(1) 



where p{r) is the gas density profile, T(r) is the gas 
temperature profile, G is the gravitational constant, k 
is Boltzmann's constant, m v is the proton mass, and 
fi is the mean molecular weight which we have taken 
to be 0.6. As is apparent from eq. 0, the total mass 
enclosed within a given radius is dependent on the 
local temperature and density gradients at that ra- 
dius. For this reason, the mass profiles presented in 
the following sections may, on occasion, decrease with 
radius. The "true" gas density and temperature pro- 
files can be obtained directly from the simulated data. 
Observationally, this information is often not supplied 
directly by the data. However, both the gas density 
and temperature profiles can be obtained through use 
of the standard X-ray deprojection technique. 

The X-ray deprojection technique, which is de- 
scribed in more detail by Arnaud (1988), uses the 
assumption of spherical symmetry to obtain the X- 
ray volume emissivity from the observed X-ray sur- 
face brightness. It then uses the assumptions of 
quasi-hydrostatic equilibrium within spherical shells 
in addition to virial equilibrium and isothermality 
of the dark matter distribution to obtain the de- 
sired gas density and temperature profiles. We have 
modified the deprojection technique to allow for non- 
isothermality of the dark matter distribution by in- 
cluding a velocity dispersion profile. Although the ve- 
locity dispersion profile is easy to obtain from the nu- 
merical simulations, it is very difficult to come by ob- 
servationally. Since we know that our simulated clus- 
ters have non-isothermal dark matter distributions, 
we assume knowledge of the true velocity dispersion 
profile in all the following analysis. Neglecting this 
correction is equivalent to underestimating the veloc- 
ity dispersion and can result in a systematic under 



estimation of the mass in the central regions of even 
relaxed clusters (see §6.1.). 

Owing to the strong dependence of thermal bremsstrahlung 
and line emission on the gas density (oc n 2 ), the den- 
sity profile is well constrained by the X-ray emission. 
The temperature dependence is considerably weaker 
(~ T 1 / 2 ), and consequently, the temperature profiles 
are considerably less well constrained. Examination 
of eq. [l] reveals how unfortunate this is since the mass 
is directly proportional to the temperature. The tem- 
perature uncertainties are less important with regards 
to the temperature gradient which is usually small 
compared to the gas density gradient. Fortunately, 
the gas temperatures can also be obtained through 
spectroscopic methods. Unfortunately, high spatial 
resolution has not coincided with high spectral reso- 
lution. Therefore, at least until very recently, cluster 
temperatures have been emission- weighted mean tem- 
peratures forcing the assumption of isothermality. As 
mentioned before, ASCA and AXAF should remedy 
this situation. 

With this in mind, we now estimate the masses of 
the simulated clusters using a variety of methods. As 
a baseline for judging uncertainties in the mass esti- 
mates, we begin by estimating the mass of our ini- 
tial, spherically symmetric, hydrostatic primary clus- 
ter. We then use the known gas density and tempera- 
ture profiles and eq. [l] to estimate the cluster masses 
in the evolutionary sample (§2.1.). Since we are using 
the true profiles and there are no projection effects, 
this analysis represents a best case scenario in which 
only the assumptions of hydrostatic equilibrium and 
spherical symmetry are being tested. We then analyze 



the deprojection sample (§2.2.) within the constraints 
of true cluster observations. That is, we deproject 
the synthetic X-ray images obtaining gas density and 
temperature profiles which are then applied to eq. El 
We also analyze the deprojection sample substituting 
a mean spectroscopically derived isothermal temper- 
ature for the deprojected temperature profile. Then, 
having determined the cluster gravitational and gas 
masses via deprojection, we examine uncertainties in 
the derived baryon fraction. Finally, we examine the 
uncertainties in virial mass estimates. 

6.1. The Hydrostatic Spherical Cluster 

Figure [li] shows the true mass profile for the pri- 
mary cluster in the 4:1 mass ratio merger. Also in- 
cluded in Figure |l4] are the mass profiles derived from 
applying eq. [l]to: 1) the true gas density and temper- 
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ature profiles and 2) the deprojected profiles out to a 
radius of 2 Mpc. The mass profile derived from the 
true profiles is quite accurate. The RMS variations 
for r >0.25 Mpc are only ~2%. The mass profile 
derived from the deprojected profiles are systemati- 
cally low by ~6% with RMS fluctuations of ~3%. In 
general, underestimates of the mass result from un- 
derestimates of the temperature. As discussed below, 
systematic errors in the deprojected temperature pro- 
file are difficult to avoid. 

The X-ray deprojection method requires as inputs, 
in addition to the X-ray surface brightness profile, a 
velocity dispersion, a core radius and an initial guess 
of the temperature (or pressure) at the outer edge 
of the profile. The velocity dispersion is obtained 
from spectroscopic redshifts of member galaxies. The 
core radius, r c , results from fitting a j3 model (Sarazin 
1986), 



Sx(r) = A 



1+ I - 



.-. -3/3+0.5 



(2) 



to the X-ray surface brightness profile, Sx- The 
temperature is often obtained by X-ray spectroscopic 
methods (e.g., David et al. 1993). All three will af- 
fect the resulting gas density and temperature pro- 
files, and ultimately the mass profile to varying de- 
grees. In terms of the mass profile, the estimate of 
the outer temperature is not significant as long as the 
profile extends beyond the bulk of the cluster mass. 
This is equivalent to saying that the choice of the 
outer temperature does not strongly influence the in- 
ner temperatures since it is an iterative solution. Of 
course, if the observed X-ray surface brightness pro- 
file extends <1 Mpc, as is often the case, then the 
choice of the outer temperature must be a good one. 
At 2 Mpc, the deprojection-based mass estimate is 
not sensitive to the choice of the outer temperature 
at least for these relatively concentrated clusters. 

The choice of velocity dispersion and core radius 
are more important because they define the gravita- 
tional potential within which the gas is supposedly in 
hydrostatic equilibrium. Figure |l5| a shows how errors 
in both the velocity dispersion and the core radius 
are reflected in the deprojected temperature profile. 
A 20% error in the core radius (only 1 zone in the 
simulation), results in a nearly 30% error in the cen- 
tral temperature. Fortunately, the problem is well 
localized in the core of the cluster and the error in 
the mass estimate, at least in this example, is only 



~10%. This is not to say that the mass error will not 
be large when the uncertainty in the core radius is 
large. At this time, it is important to note that the 
core radius one wishes to use in the deprojection is 
the dark matter (or gravitational) core radius. How- 
ever, we can only measure the X-ray core which, as 
shown by these simulations (§[?]) and various gravita- 
tional lensing experiments (Wu 1994; Miralde-Escude 
& Babul 1995), may be two or more times larger than 
the dark matter core. 

The temperature is even more sensitive to the ve- 
locity dispersion. A 20% error in the velocity dis- 
persion can cause a nearly 50% error in the central 
temperature, and the large errors are not localized 
to the core region as before. The corresponding er- 
rors in the mass estimates approach 40%. Errors in 
the velocity dispersion can easily exceed the conser- 
vative numbers used in this example. Typical uncer- 
tainties are ~10-25% (e.g., Girardi et al. 1993), but 
these are formal errors and do include systematic er- 
rors introduced by anisotropy in the velocity distribu- 
tions. Mergers produce bulk streaming and sustained 
anisotropic velocity dispersions both of which result in 
contamination of the velocity distribution and poten- 
tial projection effects. There is also a potential bias 
between galaxy and dark matter velocity dispersions 
resulting from dynamical friction effects. The veloc- 
ity dispersion required by the deprojection is actually 
the dark matter velocity dispersion which, according 
to numerical simulations (Frenk et al. 1995), may ex- 
ceed the observable galaxy velocity dispersion by as 
much as 20-30%. 

Figure |l5| b shows that uncertainties in the core ra- 
dius and velocity dispersion do not have a significant 
effect on the deprojected gas density profiles. 

6.2. Hydrostatic Mass Estimates: True Pro- 
files 

In this section, we estimate the cluster masses of 
the evolutionary sample assuming hydrostatic equilib- 
rium and spherical symmetry. We have azimuthally 
averaged the simulated 3-dimensional density and 
temperature data in spherical shells of <ir=50 kpc. 
The resulting profiles have been applied to eq. ^, and 
the results can be seen in Figure UQ. 



Figure 16 shows the true mass profile (dashed) 
overlaid with the mass profile derived from eq. [I] 
(solid) out to a radius of 2 Mpc. We find localized 
regions where the mass can be overestimated by more 
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than a factor of 2. These regions are associated with 
shocks in the ICM. As time passes, the shocks prop- 
agate away from the cluster core and dissipate. The 
resulting discrepancy between the true and estimated 
mass profiles likewise decreases and appears at larger 
radii. Whereas the mass is overestimated in the re- 
gion of the shock, the mass is underestimated in the 
region behind the shock by nearly 40%. The greatest 
errors in the mass determination occur shortly before 
and within 2 Gyrs after core passage. This point is 
demonstrated graphically in Figure [17| 

Figure ^ shows the evolution of the mean abso- 
lute error in the hydrostatic equilibrium-based mass 
estimate. We have averaged the absolute fractional 
discrepancy between the estimated and true mass pro- 
files for r >0.25 Mpc. Also included is the RMS of 
the errors over the same range in radius. Peak-to- 
peak errors are approximately three times the RMS 
errors. Again, the errors are most significant, 25% 
mean and 40% RMS, in the first 2 Gyrs after core 
passage. This is not surprising since it is during the 
first 2 Gyrs that the gas dynamics play the greatest 
role and when the conditions for hydrostatic equilib- 
rium are most severely violated. It is also not surpris- 
ing then that the shape of Figure [l?] mimics that of 
Figure [n] which shows the evolution of the mean gas 
dynamics for each of the simulations. At later epochs, 
the shocks and bulk flows have somewhat dissipated 
and both the mean and RMS errors are of order 10%. 
There is a slight trend toward systematically low mass 
estimates during the later stages of the merger. 

6.3. Hydrostatic Mass Estimates: Deprojected 
Profiles 

In this section, we estimate the cluster masses of 
the deprojection sample based on the deprojected gas 



density and temperature profiles. As stated in §5.1 



the deprojection requires an estimate of the tempera- 
ture at the outer edge of the profile, the velocity dis- 
persion of the dark matter, and the core radius. The 
outer temperature was chosen to be the true tempera- 
ture obtained from the simulation data. The velocity 
dispersion is a deprojected 1-dimensional dispersion 
profile along the "observer's" LOS at the appropri- 
ate projection angle. The core radius was obtained 
by fitting a (5 model (e.g., Sarazin 1986) to the az- 
imuthally averaged synthetic X-ray surface brightness 
at the given projection. A representative velocity dis- 
persion and the fitted core radii can be seen in Table 





In order to more closely mimic the observational 
analysis, we have chosen X-ray morphologies that ap- 
pear nominally relaxed (see Figure |2|). In those im- 
ages containing obvious substructure such as multiple 
peaks, we subtracted the secondary peaks from the 
image before obtaining the X-ray surface brightness 
profile. This is similar to the analysis done by Henry 
et al. (1993) on A2256 in which the region of the X- 
ray image containing the subcluster peak was simply 
avoided. 

The results of the deprojection can be seen in Fig- 
18 and [l9]. In Figure [l^, we display the true 



and deprojected gas density profiles for each epoch. 
As one can see, the gas density is well constrained 
by the X-ray surface brightness. Errors are typically 
small and localized in the core where resolution ef- 
fects (i.e., smoothing of the profiles) are most signif- 
icant. In Figure |l9|, the discrepancies between the 
two temperature profiles can often be large particu- 
larly in the core where they can approach a factor 
of two or more. These errors result from both the 
non-equilibrium conditions and the systematic errors 



of the deprojection discussed in §6.1 

Figure ^ shows the mass profiles derived from the 
deprojected temperature and density profiles overlaid 
with the true mass profiles. As expected, the errors 
are considerably larger than those obtained from the 
true profiles (§3.2.). Peak-to-peak errors approach a 
factor of 3. 

In the early merger examples shown here, the mass 
estimates are systematically high by 50% and >100% 
at —0.25 and 0.0 Gyr, respectively. These errors are 
largely due to an overestimation of the velocity dis- 
persion which is contaminated by the infall velocity 
of the subcluster. In the examples presented here, 
the effect is enhanced because the merger axis is pro- 
jected out of the plane of the sky. In Figure ptj[b, the 
projection is greater than in Figure p0|a, and the error 
in the mass determination is correspondingly larger. 

In the examples depicting the period shortly after 
core passage (Figure p0| c and d), there do not appear 
to be any strong systematic errors in the mass esti- 
mates. This may result because two competing effects 
are cancelling out. In both Figures |20| c and p0|d, the 
X-ray core radii have increased relative to the dark 
matter core. This should result in a underestimation 
of the central temperature and mass. However, this 
effect is compensated for by the overestimation of the 
velocity dispersion which raises the temperature and 
mass estimate. Projection effects can be severe dur- 
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ing this period. It is difficult to generalize, but merger 
remnants in the plane of the sky will tend to have their 
masses underestimated by the standard deprojection 
method while those significantly out of the plane will 
tend to have their masses overestimated. 

For t >2 Gyrs, the examples presented here tend 
to have their masses systematically underestimated. 
In both Figures ppj e and pp[ f, the X-ray core radii are 
significantly larger than the dark matter core radii 
(nearly a factor of 2 in pO| f) causing the central tem- 
peratures and masses to be severely underestimated. 
Because of the prolate structure of the merger rem- 
nant, the X-ray cores are subject to projection ef- 
fects such that they will appear larger when viewed 
perpendicular to the merger axis. Projection effects 
in the velocity dispersion tend to be minimized dur- 
ing the later merger stages. There is some residual 
anisotropy in the velocity dispersion (^5-15%) even 
after 5 Gyrs, but it is small compared to the velocity 
anisotropy early in the merger (~25-50%). 

Figure |l| shows the evolution of the mean emissivity- 
weighted error in the deprojected mass estimate. We 
have averaged the fractional discrepancy between the 
deprojected and true mass profiles for r >0.25 Mpc 
and weighted it by the X-ray surface brightness at 
that radius thereby accentuating uncertainties in the 
brighter central regions of the clusters where one is 
more likely to make observational mass estimates. Of- 
ten, we find post-merger mass estimates systemati- 
cally low by ^50%. As a comparison, we have also 
included the mean absolute error at the same epochs. 
The mean is taken over the same range in radius, but 
is not emission-weighted. Comparison with Figure 
[l7| shows that the mass errors are not so exclusively 
linked to the gas dynamics as in §6.2. This would 



seem to imply that significant systematic effects are 
being introduced by the deprojection technique. The 
mean emissivity-weighted errors show a trend toward 
significant underestimation of the mass even after 5 
Gyrs. The discrepancy between the mean absolute 
errors and the mean emissivity-weighted errors shows 
that errors are typically larger at small radii. As dis- 
cussed above, the inflated X-ray cores cause a system- 
atic underestimation of the central temperature and 
an underestimation of the mass. 

6.4. Isothermal Spectroscopic Temperature 

Since the errors in the mass estimates are tied so 
closely to the gas temperatures, it makes sense to 
use spectroscopically determined temperatures rather 



than the deprojection based temperatures whenever 
possible. However, even the spectroscopically deter- 
mined temperatures are poorly determined (e.g., David 
et al. 1993), and the limited spatial resolution requires 
the assumption of isothermality. Since merging clus- 
ters are not isothermal (§[!]) we need to assess the un- 
certainties introduced by that assumption. We have 
therefore recalculated the mass estimates as before, 
using the deprojected density profiles, although this 
time we have replaced the deprojected temperature 
profiles with mean emissivity-weighted temperatures. 
This method helps to reduce some of the systematic 
errors noted previously. However, since the density 
and temperature profiles are no longer self-consistent, 
considerable noise is added to the derived mass pro- 
files. 

In Figure |22|, we have plotted the same quantities 
as in Figure |2l| only now for the isothermal analy- 
sis. We find that the mean emissivity-weighted errors 
have been reduced significantly by inclusion of the im- 
proved temperature estimates. The errors are always 
less than 50% and for >2 Gyrs after core passage they 
are less than 25%. We also note that the mean abso- 
lute errors are, at least for r >0.25 Mpc, somewhat 
larger than in Figure |l]. This is due to the increased 
noised induced by restricting the mass calculation to 
isothermal conditions. Peak-to-peak errors, can be 
quite large. Typically, they are of order 50-100%, but 
in extreme cases, they can exceed a factor of 3. 

6.5. A Comparison with Previous Work 

Schindler (1996) quotes an RMS uncertainty of 
15% for mass estimates within the central 2 Mpc of 
a cluster for her merger simulations. However, she 
does note several instances where the errors can be 
considerably larger. In the case of isothermal clus- 
ters, deviations are between 20% and 35%. Near the 
radius of shocks, cluster masses can be overestimated 
by more than 100%, and substructure can cause se- 
vere underestimations at the radius of the substruc- 
ture. She also notes that in some instances errors in 
the mass estimates can be greater at radii <0.5 Mpc 
(20-30%) than at larger radii (10-20%). 

We find that these results are in general agreement 
with our own although our uncertainties tend to be 
somewhat larger, and we find an epochal dependence 
in the errors. Our quoted errors tend to be larger 
because we used deprojected density and tempera- 
ture profiles (§ |6.3.| ) or mean isothermal temperatures 
(§6.4.) whereas she used the deprojected density pro- 
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file along with the true temperature profile. Also, we 
have weighted our errors by the X-ray surface bright- 
ness at the given radius thus emphasizing the more 
observationally important inner regions of the cluster. 
Most clusters do not have X-ray emission detected 
beyond 0.5-1.0 Mpc let alone 2-3 Mpc. Finally, we 
find the greatest uncertainty in mass estimates made 
within the first 2 Gyrs after core passage when equi- 
librium conditions are most severely violated. 

6.6. Implications for the Baryon Fraction 

Uncertainty in the baryon mass fraction has two 
sources. They are the uncertainty in the total mass 
estimate, as described above, and the uncertainty in 
the total gas mass. The total gas mass can be ob- 
tained simply by integrating the deprojected gas den- 
sity profiles. The similarity between the deprojected 



density profiles and the true profiles (see Figure 18), 
would seem to indicate that the total gas mass is well 
constrained by the X-ray emission. We find errors in 
the integrated gas mass to be typically of order 5%. 
There are a few exceptions when errors exceed 10% 
at a given radius. 

Since the uncertainties in the total gas masses are 
small relative to those in the total mass, the primary 
source of uncertainty in the baryon fraction resides 
in the total mass estimate. Figure ^3| shows the ratio 
of the observed (i.e., deprojection-based) baryon frac- 
tion and the true baryon fraction (solid line) within a 
given radius for the deprojection sample. As a com- 
parison of the relative uncertainties, we have also in- 
cluded the ratio of the integrated deprojection-based 
gas mass to the true gas mass (dotted line). In both 
ratio of unity indicates an agreement between 
"observations" and "truth." 

In the examples of the early merger stages (Fig- 
ure ^3| a and b), the baryon fraction is systemati- 
cally underestimated by ^50%. Recall that the to- 
tal mass was overestimated. In the examples of the 
later merger states (Figure f23| c, e, and f), there is a 
systematic overestimation of the baryon fraction, par- 
ticularly in the central regions (r <1.0 Mpc). The de- 
gree of overestimation, often greater than 50%, could 
help to reduce the discrepancy between observed and 
theoretically derived baryon fractions. As an exam- 
ple, White & Fabian (1995) find the baryon fractions 
in a sample of 13 clusters to range between 0.1-0.22 
within 1 Mpc while Big Bang nucleosynthesis calcu- 
lations predict a value of ~0.06 for an f2=l universe 
(Olive et al. 1990). Both the observed and theoreti- 



cal values assume H o =50 km s _1 Mpc -1 . It should be 
kept in mind that much of the discrepancy observed 
in Figure ^3| is due to the systematics of the deprojec- 
tion technique as revealed in the mass estimates. If we 
calculate the baryon fraction as a function of radius 
using the deprojected gas density and the total mass 
estimates derived as in § |6.2j (i.e., using true density 
and temperature profiles), then the RMS uncertain- 
ties in the baryon fraction are typically less than 15% 
except within the first 1 Gyr after core passage in 
which case they can exceed 40%. 

6.7. Virial Mass Estimates 

As an illustrative example of the relative errors in- 
volved in optical velocity-based mass estimates, we 
have calculated the virial masses of each cluster sys- 
tem in the deprojection sample. The virial mass is 
given by 

M vT=^<(iy 1 (3) 

where a v is the LOS velocity dispersion, <l/r> _1 
is the harmonic mean radius (e.g., Oegerle & Hill 
1994), and G is the gravitational constant. It is im- 
portant to remember that Mvt is the total cluster 
mass whereas the previous X-ray analysis provided 
the mass enclosed within a given radius (eq. |l|). In 
this analysis, we calculated the velocity dispersion 
and harmonic mean radius for 100 draws of 100 ran- 
domly selected particles within 2.5 Mpc of the pro- 
jected center of mass. Applying this analysis to the 
initial, virializcd primary cluster (4:1 mass ratio), we 
found MyT=3.9±.5 x 10 14 M Q compared to a true 
mass of 4.0 x 10 14 M Q for a 2.5% error. The results of 
the same analysis applied to the deprojection sample 
can be seen in Table § 

From Table it is apparent that the errors are 
greatest when the velocity dispersion is most strongly 
contaminated by the infall velocity. This occurs early 
in the mergers (near core passage), when the merger 
axis is projected out of the plane of the sky, and 
for lower mass ratio mergers where the probability of 
drawing a subcluster particle is higher. The minimum 
errors occur for high mass ratio mergers, mergers in 
the plane of the sky, and in the later stages of the 
merger when the system is returning to equilibrium. 
This result differs from that of Bird (1995) who finds 
that correcting for substructure in virial mass esti- 
mates more strongly affects the mean galaxy separa- 
tion rather than the velocity dispersion. 
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Obscrvationally, there are more sources of error 
than are addressed in this analysis. We do not in- 
clude contamination by foreground and background 
galaxies, a galaxy/dark matter velocity bias, nor do 
we include measurement errors in the individual par- 
ticle velocities. It should also be noted that we do not 
attempt to remove obvious substructure which could 
somewhat reduce the errors. Identification of sub- 
structure is strongly dependent on a complete sam- 
pling of both the galaxy spatial and velocity distribu- 
tions and is therefore very difficult in most clusters. 

7. SUMMARY 

In this paper, we have outlined the many observ- 
able consequences of cluster-subcluster mergers find- 
ing evidence of merger activity in a variety of clus- 
ter observations. We have also investigated the ef- 
fects of recent merger activity and the resulting non- 
equilibrium conditions on the accuracy of equilibrium- 
based cluster mass estimates and estimates of the 
cluster baryon fraction. 

In the X-ray emission, we not only find evidence 
of obvious substructure (i.e., multiple peaks), but we 
also find more subtle evidence of non-equilibrium con- 
ditions such as isophotal twisting, ccntroid shifting 
and the enlargement of the X-ray core relative to the 
dark matter distribution. Elongation of the X-ray 
morphology is also a consequence of mergers although 
it is relatively long-lived and therefore may not be a 
sensitive indicator of recent merger activity. 

The X-ray emission also contains valuable temper- 
ature information. We find that one of the strongest 
indicators of recent merger activity is in the temper- 
ature morphology of the X-ray emitting gas. Mergers 
tend to result in rather large temperature inhomo- 
geneities (>few keV) on scales <0.5 Mpc. This result 
is quite consistent with recent high resolution temper- 
ature maps of several rich clusters (e.g., Briel & Henry 
1994; Henry & Briel 1995, Markevitch et al. 1994; Ar- 
naud et al. 1994), all of which have other signatures 
of recent merger activity. 

The gas dynamics also result in several observable 
consequences. We have presented a LOS velocity map 
as an example of what might become available with 
the launch of ASTRO-E. The merger results in a very 
complex velocity structure that persists for several 
Gyrs after the merger. The specifics of the observa- 
tional analysis need to be addressed in more detail in 
future work. 



We have added quantitative support to the idea 
that bulk gas flows resulting from mergers can pro- 
duce the wide-angle tailed (WAT) radio source mor- 
phology. Although not specifically addressed in this 
work, we also believe that bulk flows may result in 
the disruption of cooling flows and that this is true 
even in high mass ratio mergers. McGlynn & Fabian 
(1984) proposed that mergers of nearly equal mass 
clusters might disrupt cooling flows. We also expect 
mergers to disrupt the Ha filaments which are often 
associated with cooling flows (sec Baum 1992) and 
thus may explain the results of Donahue, Stocke, & 
Gioia (1992) which show a reduction in the number of 
cooling flows (as indicated by Ha filaments) in recent 
epochs. We also link mergers, via shocks, turbulence 
and bulk flows, to the formation of Mpc-scale radio 
halos. We further contend, based on the decay time 
for the ICM dynamics, that halos may be relatively 
long-lived compared to the canonical synchrotron life- 
time of 10 8 years. 

Mergers may also influence the evolution of indi- 
vidual galaxies. As described here, mergers generate 
expansive shocks in the ICM through which individ- 
ual galaxies may pass. In doing so, they will experi- 
ence an increase in ram pressure which can induce a 
burst of star formation changing the observed galaxy 
colors and stellar content. The ram pressure may also 
strip the galaxy of its gas content reducing the po- 
tential for future star formation and thus altering the 
galaxy's observed color and morphology. 

The dynamics induced by the merger were shown 
to affect the equilibrium-based cluster mass estimates. 
Even in a best case scenario, (i.e., knowledge of the 
true temperature and density profiles), RMS errors 
in the mass estimates can approach 40% in the first 2 
Gyrs after core passage because of shocks and bulk 
flows in the ICM. At later times, RMS errors are 
reduced to ^10%. Local errors in the mass esti- 
mates, particularly in the early stages can be in ex- 
cess of 100%. Using density and temperature pro- 
files obtained from the standard X-ray deprojection 
technique adds potentially severe systematic errors. 
These can be somewhat reduced by using a spectro- 
scopically derived temperature rather than the de- 
projected temperature profile. Still, the discrepan- 
cies in the mass estimates can be large, approach- 
ing 50% within 2 Gyrs of core passage and 25% for 
later epochs. These results are consistent with those 
of Evrard (1994) and somewhat larger than those of 
Schindler (1996). 



16 



Of course, it should be recognized that we have 
only analyzed a limited sample of merger remnants 
and extrapolation to the Universe at-large is rather 
tricky. First, the uncertainties can vary markedly for 
any given cluster depending on the details of the anal- 
ysis and the radial extent of the X-ray data since there 
is a radial dependence on the fractional mass discrep- 
ancies. Also, projection effects can be severe. Second, 
the epochal dependence on the mass estimate discrep- 
ancies requires knowledge of merger rates. That is, if 
major mergers occur more often than every 2 Gyrs, 
then finding a truly relaxed cluster may be very dif- 
ficult. It has been suggested by Edge, Stewart & 
Fabian (1992) that any given cluster will experience 
a merger every 2-4 Gyrs. 

The errors in the mass estimates translate directly 
into errors in the observed baryon fraction within 
clusters of galaxies. The tendency to underestimate 
the total mass throughout the post-merger phase may 
account for some of the discrepancy between the ex- 
pected and observed baryon fractions in some clus- 
ters. 
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Figure Captions 

Fig. ^ Synthetic X-ray surface brightness images 
of a 4 :1 mass ratio merger (the evolutionary sample, 
§ 2.1. ) at six epochs. The epoch is is listed above each 
panel and is relative to the time of core passage. The 
subcluster enters from the top of the panel and pro- 
ceeds downward. Contours are uniformly spaced in 
the logarithm and constant from one panel to the 
next, with dynamic range log(S ma:E /S m i n )=5.4. Each 
panel is 6.5 Mpc on a side. 

Fig. Synthetic X-r ay s urface brightness of 

the deprojection sample (§ |2.2j ). Panels a-f corre- 
spond directly to Mergers 1-6 in Table [j]. Each 
panel is 6.5 Mpc on a side. The contours are uni- 
formly spaced in the logarithm and relative to the 
peak surface brightness in each panel. The dynamic 
range in surface brightness represented in each panel 
is bg(S ma:c /S min )=5.0, 5.6, 4.5, 4.6, 4.25, 4.55. 

Fig. [|: The relative expansion of the gas core 
(dashed) with respect to the dark matter core (solid) 
as a result of heating and gas motions generated dur- 
ing the 4:1 mass ratio merger, a) Initially, the gas 
and dark matter cores are of comparable size, b) At 5 
Gyrs after core passage, the gas core has nearly dou- 
bled in size where the dark matter core is virtually 
unchanged. 

Fig. ^ The result of fitting elliptical isophotes to 
the X-ray surface brightness image in Fig. ||d. The di- 
mensions of the box are 3 Mpc on a side. The merger 
axis runs vertically. This figure shows the elongation 
of the X-ray distribution along the merger axis, a ra- 
dially dependent ellipticity, and a significant centroid 
shift of - 150 kpc. 

Fig. ^ The emissivity-weighted temperature dis- 
tribution for the evolutionary sample (§ [2.1.| ). Solid 
contours are labeled with the temperature in keV. 
The X-ray surface brightness is shown as dotted con- 
tours that are uniformly spaced in the logarithm over 
the same dynamic range as in Figure [j]. Each panel 
is 4 Mpc on a side. 

Fig. ^ The emissivity-weighted tem perat ure dis- 
tribution for the deprojection sample (§ [2.2.| ). Solid 
contours are labeled with the temperature in keV. 
The X-ray surface brightness is shown as dotted con- 
tours that are uniformly spaced in the logarithm over 
the same dynamic range as in Figure ^. Panels a-f 
correspond directly to Mergers 1-6 in Table [|. Each 
panel is 4 Mpc on a side. 



Fig. |?[ A slice in temperature (filled contours) 
along the merger axis through the cores of a 2:1 mass 
ratio merger. The merger axis runs vertically with 
the slightly cooler subcluster at the top and moving 
downward. The epoch is 0.25 Gyrs before core pas- 
sage. The contours are evenly space between 0.5 and 
8.5 keV. Black designates low temperatures. White 
designates high temperatures. The image dimensions 
are 6.5 x 11.2 Mpc. 

Fig. The gas dynamics, a) The observable 
emission- weighted mean LOS velocity field for Merger 
4 of the deprojection sample (§2.2.). b) The emission- 
weighted LOS velocity dispersion. In both (a) and 
(b), the velocity contours (solid) are labeled with the 
velocity in km s _1 . In (a), the velocities are in the 
rest frame of the primary cluster's dark matter core, 
and negative velocities are approaching the observer. 
The dotted contours are the X-ray surface brightness. 
Both panels are 2 Mpc on a side. 

Fig. ^: The evolution of the peak gas velocity 
within a 2 Mpc box centered on the dark matter den- 
sity peak for mergers having mass ratios of 1:1 (+), 
2:1 (*), 4:1 (o), and 8:1 (A). 
Fig. 
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The evolution of the width of bulk flows 
through the cluster core for mergers having mass ra- 
tios of 1:1 (+), 2:1 (*), 4:1 (o), and 8:1 (A). The flow 
width is defined as the width at the primary dark 
matter core at which the flow velocity drops to 0.5 of 
the peak or at which the flow direction varies from 
that of the merger axis by more than 25°. 
Fig. 
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The evolution of the mean gas velocity 
within a 2 Mpc box centered on the dark matter den- 
sity peak for mergers having mass ratios of 1:1 (+), 
2:1 (*), 4:1 (o), and 8:1 (A). 
Fig. 
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The evolution of the median ram pres- 
sure experienced by the N-body particles during the 
8:1 mass ratio merger. Subcluster particles (o), on 
average, experience an impulsive increase in ram pres- 
sure near the time of core passage. A smaller but still 
sizable fraction of primary cluster particles (A) also 
experience a large increase in ram pressure as they 
pass through shocks generated during the merger. 
Fig. 13: Mach number contours at 0.5 Gyrs after 



core passage in the 2:1 mass ratio merger. Shocks are 
identified as the two opposing C-shaped structures 
(closely spaced contours). The cluster core is located 
between these structures. The region depicted is 6.5 
x 11.2 Mpc. 
Fig. 



14: Mass Estimates of the initial hydrostatic 
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cluster. The solid line is the true mass profile. The 
dotted line represents the mass profile derived using 
Eqn. |l| and the true gas density and temperature 
profiles. The dashed line represents the mass profile 
derived using Eqn. |l| and the deprojection-based gas 
density and temperature profiles. 
Fig. 
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A comparison of the systematic errors 
associated with the X-ray deprojection of an isother- 
mal cluster in hydrostatic equilibrium, a) Errors in 
the temperature profile induced by errors in the core 
radius (o) and velocity dispersion (o). Filled symbols 
represent 20% overestimates of the quantity. Open 
symbols represent 20% underestimates of the quan- 
tity. The solid line represents the profile derived us- 
ing the best estimates of the core radius and velocity 
dispersion. The "true" gas temperature is 4.6 keV. b) 
The corresponding errors induced in the gas density 
profile. 
Fig. 
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A comparison of the true (dashed) and 
equilibrium-base mass profiles (solid) for the evolu- 
tionary sample (§2.1.). The evolutionary sample is a 
4:1 mass ratio merger at the six epochs listed above 
each panel. All times are with respect to core pas- 
sage. The equilibrium-based mass profile was derived 
using Eqn. |l] and the true (azimuthally averaged) gas 
temperature and density profiles. 

Fig. |l7]: The mean absolute (o) and RMS errors 
(A) in the equilibrium-based mass estimates of the 
evolutionary sample (§2^L, Fig. O) as a function of 



time relative to core passage. The fractional error is 
averaged over radii from 0.25 to 2.0 Mpc. 

Fig. 18: A comparison of the true (solid) and 



deprojected (dotted) gas density profiles for the de- 
projection sample (§2.2.). Panels a-f correspond di- 
rectly to Mergers 1-6 in Table [|. All profiles have 
been normalized. 

Fig. |l9| : A comparison of the true (solid) and 
deprojected (dotted) g as te mperature profiles for the 
deprojection sample (§ |2.2 j ) . Panels a-f correspond 
directly to Mergers 1-6 in Table [j]. 

Fig. 20: A comparison of the true (solid) and 



deprojection-based ( dott ed) mass profiles for the de- 
projection sample (§2.2.). Panels a-f correspond di- 
rectly to Mergers 1-6 in Table |. 

The mean emissivity-weighted error (o) 



Fig. 21 



and mean absolute error (A) for the deprojection- 
based mass estimates (Fig. ^0|) as a function of time 
relative to core passage. Since various merger pa- 
rameters are represented in this sample, this figure 



cannot be viewed as a true time sequence although 
the magnitude of the uncertainties depicted here are 
considered to be representative of the various epochs. 
The fractional errors have been averaged between 0.25 
and 2.0 Mpc for each profile in Fig. |2(]. 

Fig. 22: The mean emissivity-weighted error (o) 
and mean absolute error (A) for the deprojection- 
based mass estimates in which the deprojected tem- 
perature profile has been replaced by the mean emissivity- 
weighted temperature. As in Fig. [2l], this is not a 
true time sequence, and the fraction error is averaged 
between 0.25 and 2.0 Mpc for each mass profile. 
Fig. 
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The error in the cluster baryon fraction 
based on the deprojected gas density profiles and the 
equilibrium/deprojection-based mass profiles. The 
dotted line represents the uncertainty in the baryon 
fraction assuming deprojection-based gas mass esti- 
mates and knowledge of the true mass profile. Unity 
implies a good agreement (dashed line). The solid 
line represents the uncertainty in the baryon fraction 
estimate assuming both total masses and gas masses 
based on deprojection. Panels a-f correspond directly 
to Mergers 1-6 in Table fjj 
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Initial Cluster Parameters: 


Deprojection Sample 






ID 


Mass Ratio 


Primary Mass 


Epoch 


P.A. 


a v 


<T> 


V 






(10 14 M ) 


(Gyrs) 


Degrees 


(km s- 1 ) 


(KeV) 


(kpc) 


1 


2:1 


2.0 


-0.25 


45 


790 


2.8 


245 


2 


4:1 


4.0 


0.0 


90 


1300 


8.1 


162 


3 


8:1 


8.0 


0.5 


45 


1055 


9.5 


411 


4 


2:1 


2.0 


1.0 


45 


590 


2.5 


142 


5 


4:1 


4.0 


2.0 





660 


4.5 


333 


6 


8:1 


8.0 


5.0 





995 


9.6 


492 



Table 1: Column 1 is the merger ID. Column 2 is the cluster mass ratio. Column 3 is the primary cluster mass. 
Column 4 is the merger epoch relative to the time of core passage. Column 5 is the projection angle of the merger 
axis with respect to the plane of the sky. Column 6 is the LOS velocity dispersion. Column 7 is the mean emissivity- 
weighted temperature within 1 Mpc of the X-ray maximum. Column 8 is core radius obtained from fitting a (3 
model (eq. 



Virial Mass Comparison 
ID True Mass M VT 



(10 14 M ) (10 14 M Q ) 



1 


3.0 


5.1±0.6 


2 


5.0 


10.0±0.15 


3 


9.0 


10.2±0.15 


4 


3.0 


4.0±0.7 


5 


4.7 


4.4±0.7 


6 


8.6 


8.3±1.1 



Table 2: Column 1 is the merger ID within the deprojection sample. Column 2 is the true total mass. Column 3 
is the virial mass estimate. 
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Fig. 3.- 
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Fig. 8.- 
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Fig. 9.— 
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Fig. 10.— 
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Fig. 13.— 
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Fig. 17.— 
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Fig. 21.— 
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Fig. 22.— 
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